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■ The mechanisms of stability or instability in the strained alloy film growth are of intense current 
interest to both theorists and experimentalists. We consider dislocation- free, coherent, growing 
alloy films which could exhibit a morphological instability without nucleation. We investigate such 
strained films by developing a nonequilibrium, continuum model and by performing a linear stability 
analysis. The couplings of film-substrate misfit strain, compositional stress, deposition rate, and 
growth temperature determine the stability of film morphology as well as the surface spinodal 
decomposition. We consider some realistic factors of epitaxial growth, in particular the composition 
dependence of elastic moduli and the coupling between top surface and underlying bulk of the film. 
The interplay of these factors leads to new stability results. In addition to the stability diagrams 
both above and below the coherent spinodal temperature, we also calculate the kinetic critical 
thickness for the onset of instability as well as its scaling behavior with respect to misfit strain and 
deposition rate. We apply our results to some real growth systems and discuss the implications 
related to some recent experimental observations. 
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I. INTRODUCTION 



The techniques and technology related to fabricating large single crystals of elemental materials like silicon (Si) or 
compounds like gallium arsenide (GaAs) has been perfected for decades now. Polished wafers of such single crystalline 
semiconductor materials are now routinely used as substrates, on which a series of thin layers are grown to create a 
desired heterostructure. Thin films of increasing complexity are now being fabricated due to constraints on operating 
conditions and performance of electronic and optical devices. 

When the growth of the superimposed material layers is based on the atomic pattern of the substrate underneath, 
as if it is the extension of substrate structure, then the process is referred to as an epitaxial growth. Heteroepitaxy 
is the technique in which a thin film of a different material is grown on a substrate (as for example in a layer 
growth of AlAs on GaAs), and in which the deposited layer has the same crystal orientation as the substrate. Many 
materials (especially III-V compounds), though composed of different atoms, have identical crystal structures and 
' nearly identical lattice spacing. A match of lattice spacing is important in epitaxial growth since it can minimize and 
possibly eliminate local strains in the growing film, and thus produce good quality thin solid layers. For example, one 
can in principle lattice-match an InGaAsP quaternary alloy of appropriate composition to either InP or GaAs, and 
tailor-make a semiconductor film which will have the electronic band gap within the range between 0.8 eV and 1.7 
eV. 

Although the alloy composition can be tuned to obtain a desirable band gap, the resulting device is useful only if the 
alloy composition remains homogeneous. In a typical epitaxial growth technique like molecular beam epitaxy (MBE), 
the stability of the growing film can depend on many variables. The homogeneity of the growing film clearly depends 
on the surface mobility of the material being deposited, the substrate temperature T, and the deposition rate v. 
Further, different sized atoms that are being deposited to fabricate an alloy film induce stresses in the substrate which 
in turn affect the behavior on the surface. In fact, the process of epitaxial growth is a nonequilibrium phenomenon 
in which surface, interfacial and elastic energies along with surface diffusion and alloy segregation play active roles. 
During the nonequilibrium growth of alloy films, elastic stress may be relieved through nucleation of dislocations and 
other defects. Inhomogeneities can also arise on account of two potential instabilities, from which the generation 
of dislocations can be avoided and one can obtain high quality thin film materials. These are the morphological 
instability and alloy segregakkm instability. 

Morphological instabilitytiHu can be understood through two features: A solid film adsorbate on a solid substrate 
creates a uniaxial stress, and the driving force for surface diffusion of matter is proportional to the local strain energy 
density. If a small corrugation occurs on the surface of a uniaxially strained solid, local strain energy concentration is 
created at its valleys. As a result diffusion of matter occurs along the surface from valleys to peaks of the corrugation. 
Valleys then deepen due to loss of mass which leads to a further increase of local strain energy driving the instability. 
Surface energy plays a stabilizing role here: as corrugations deepen, surface area increases which has an associated 
energy cost inhibiting the growth of instability. n 

The alloy segregation instability is driven by alloy thermodynamics .□ It is typically modelled by a Ginzburg-Landau 
type free energy (the order parameter cf> is related to the binary alloy concentration) in which the quadratic term has 
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a coefficient proportional to (T — T c ). This term changes sign for T < T c and leads to an instability which drives a 
homogeneous mixture to segregate. This mechanism is called spinodal decomposition and is applicable for not too 
asymmetric mixtures within the classical spinodal. For asymmetric mixtures between the spinodal and coexistence 
curves, the mixture is metastable and requires large amplitude fluctuations to segregate via a nucleation mechanism. 
In this paper we consider strained film growth without nucleation, and assume the growing film to be a coher- 

Stability of such films have been considered previously both 
d experimentally. The systems explored experimentally include 
Other interesting experimental studies include: single component 
self-organised GaAs islands,c2l etc.; but these are not as relevant 
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films (Ge on a Si substrate) 

here. Many of these experimental papers show that stressed growing films can develop non-planar morphologies 
without dislocations or nucleation, for both symmetric binary alloy films (50-50 mixture) as well as asymmetric films. 
The onset of instability is measured through a kinetic critical thickness h c (film thickness at which surface roughness 
first appears), and the characteristic wavelength of the instability has also been determined by experiments. 

Early theoretical work focused on explaining nonplanar surface morphologies through the competition between 
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surface energy and elastic strain energy arising from lattice mismatch e = (a,f 
film lattice parameter and a s is the substrate lattice spacing. The important 
the film lattice parameter a* , represented by the solute expansion coefficie 
concentration field, were added in a kinetic treatment of the growing filmJ 
e and 77 leads to sensitivity of the morphological development to the sign of the misfit e, i.e. 
regions of the growing film behave differently. Another important factor, the compopitkm dependence of the elastic 
moduli, can also impact on the stability of the film and its sensitivity to the sign of e.E3o Experimental observations 
can provide a sensitive test of these theories in connection with the e-asymmetry in the stability diagram of the 
growing film, the coupling between morphological instability and alloy segregation, as well as the behaviors of critical 
thickness. 
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An additional factor is the differing atomic mi 
has been investigated in recent theoretical work 



of the two (or more) species on the growing film surface, and 
This consideration is based on recent experimental results of 
atomic mobilities on surface. E.g., thej-activation energy of step mobility for Ge on Ge(001) surface was measured to 
be smaller than that of Si on Si(001),c3 corresponding to much larger atomic mobility for Ge. However, during the 
MBE growth for semiconductors (e.g., Si-Ge system or III-V alloys), dimerized atoms are quickly formed on (001) top 
surface when new materials are deposited, leading to a typical phenomenon of surface reconstruction^ These dimers, 
which may be comprised of unlike alloy components or like atoms, diffuse_uver the surface with diffusion processes 
and properties more complicated than and different from that of monomers. Ej For example, in Ge-op^i submonolayer 
system the surface pattern with Si-Ge mixed dimers rather than Ge-Ge dimers can be observed.EJ A mixed Si-Ge 
dimer can diffuse as a unit, or transform into a pure dimer between like atoms (Si-Si) due to atomic exchange with 
underlying substrate atom and can also return back to the mixed form due to re-exchange process Ml 45 ! Therefore, 
different alloy components deposited on surface may or may not diffuse independently. In the following studies we 
consider an effective surface mobility (or effective surface diffusivity D s ) to describe the processes of surface diffusion 
and decomposition in continuum approximation. 

In this paper our starting point is the work of Leonard and Desai. It is reviewed in Ref. ^ and the basic model 
on which some refinements are made is described in Ref. ^0[ Here we further develop the model to a more general 
case, and more importantly, we consider two key features in the epitaxial growth systems: The first one is the 
coupling between top surface and underlying bulk of the growing film, which leads to a new dynamic equation for 
the evolution of surface compositional variable (4> s )] the second one is the composition dependence of all the three 
elastic constants, including the Young's modulus E, shear modulus [i and Poisson ratio v (= E/2/i — 1). In Ref. |2^, 
the shear modulus was assumed to be constant and the surface-bulk coupling was neglected in dynamical equations. 
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■can lead to stability results different from that of the previous theoretical work for strained 
B2I and more complicated and richer properties of the system can be obtained. Here we 
are interested in the growing films with depositiop^rate v ^ 0, while the results for isothermal annealing alloy films 
without growth (v = 0) are presented elsewhere.cll The details of our model for the growth of thin alloy films are 
described in Sec. II. In Sec. Ill we present the solution of mechanical equilibrium for the system (More general results 
for mechanical equilibrium solution are given in Appendix |A]). In Sec. IV we use the results of Sec. [II to derive 
self-contained dynamic equations for morphological and compositional perturbations, which we then use to perform a 
linear stability analysis; this leads to the characteristic equation for perturbation growth rate. The formulae of elastic 
free energy in Fourier space, which are needed in the derivation of Sec. IV, are given in Appendix |b|. The results 
of stability analysis for both the composition independent and composition dependent elastic moduli, including the 
stability diagrams, are presented in Sec. |v|. The behavior of kinetic critical thickness is presented in Sec. VI. Finally 



we conclude with a summary of salient features of our results in Sec. VII 
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II. MODEL 



Let us consider a strained alloy system composed of a semi-infinite substrate occupying the region z < y) and 
a Ai^x^x binary or pseudo-binary alloy film in the region y) < z < h(x, y, t). Here £(x, y) refers to the vertical 
position at which the film-substrate interface is located, and h(x, y, t) is the surface height variable of the growing film. 
Usually, CLx, y) is assumed to be constant, or for simplicity 0, corresponding to a planar film-substrate interface fixed 
at z = O.Ej The local composition within the film can be denoted by a continuous variable 0(r, t) which is proportional 
to the difference in the local concentrations of two alloy constituents. Its average value is equal to 2X — 1, with X 
the alloy composition. Here we focus on the symmetric mixture film, i.e., X = 1/2 alloy for which = 0, while in the 
substrate z < £ we have = since the compositional fluctuation is usually assumed to be absent there. 

For this system of epitaxial growth, the film is assumed to be elastically isotropic and coherent with the substrate, 
without any misfit dislocations or other defects. On the top surface of the growing film, the evaporation and re- 
condensation are negligible if the deposition of material occurs under ultra-high vacuum condition, like the growth 
through MBE, and then the system evolution should correspond to conserved dynamics. We also assume that there 
is no interdiffusion between film and substrate, and the diffusion and compositional relaxation in the bulk film can be 
neglected, since the bulk mobility of film components is miniscule compared to that on the surface in typical epitaxial 
growth. The neglect of the bulk mobility implies that the composition profile within the film reflects the time-history 
of the growing film and is dependent on the deposition rate v. As well the layers within the film are buried metastable 
layers which are frozen on account of negligible mobility in the bulk. Therefore, the bulk concentration field 0b at 
some time t is equal to the surface field S of an earlier time, i.e., <fib{x , y , z , t) = (f> s (x,y,t — to) for z < h, where 
to ~ (h — z)/v with h = vt the average surface height, and then we have 

(t>b(x,y,z,t) ~ <j> a (x,y,t = z/v). (1) 

Consequently, the conserved dynamics of morphological and compositional evolution is dominated by the surface 
diffusion and surface decomposition processes, with two time-dependent essential dynamical variables: the surface 
morphology h(x,y,t) and the concentration fieldi|*t the surface 4>{x,y,z = h(x,y,t),t) = <j) s (x,y,t). The evolution 
equation describing the time dependence of h isr°H 22 l 

ln =T ^Jh +v > (2) 

where v represents the deposition rate or growth velocity, V 2 is the surface Laplacian, g = 1 + |Vft| 2 denotes the 
determinant of the surface metric, and Th = D s N s /kBTN% is the kinetic coefficients which depends on the surface 
diffusivity D s , the Boltzmann constant ks, growth temperature T, as well as the number density of atoms per unit 
surface area N s and that per unit volume N v . For the concentration field 0, the conserved dynamics leads to 

g.rvvg-A^-*), (3) 
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where A is proportional to film deposition rate v, corresponding to the tendency of surface composition to go towa 
the homogeneous phase (of X = 1/2 and = in this paper) driven by the uniform deposition flux of material,! 
and A(z — h) is zero everywhere except on the top surface z = h, where it is 1. Here the mobility is approximately 
zero in bulk, while at surface it is = ThS^ 1 with S the effective diffusion thickness of surface layer. The total free 
energy functional T in Eqs. (||) and (||) makes the evolution of variables h and coupled with each other and consists 
of three contributions: 

T = T S +T GL + T ch (4) 

with the surface energy J-" s , the Ginzburg-Landau free energy .FgL) an d the elastic energy The addition pf 

deposition noise and thermal noise terms in the dynamical equations for h and have been considered earlier.E3 
At the level of linear stability analysis, neglecting noise does not alter the results. In Eq. (||) there are two ways 
to evaluate the surface composition fluctuations. The first oneE3 is only considering the surface state and surface 
free energy, that is, evaluating the free energy T at surface z = h and then calculating the functional differentiation 
directly with respect to surface composition field (f> s (x,y,t) to obtain its evolution behavior. In this paper, we also 
take into account the coupling between surface and underlying bulk of the growing film by (i) identifying the free 
energy T as the total energy of the system, that is, from z = -co to z = ft, and not just the one evaluated at surface 
z = h, (ii) then calculating V 2 5J r /<50 as a whole, and finally (iii) evaluating its value as well as the other terms in 
Eq. (|^) at surface z = h to obtain the dynamics of surface field S . Hiis is analogous in spirit to the procedure in 
the previous study of surface critical phenomena for spin fluctuations.EU 
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The surface energy contribution T s to the free energy 
drumhead model without pinning term: 



plays a stabilizing role and can be represented by a 



(5) 

where 7 is the surface tension at the top surface assumed to be isotropic and composition independent .El The alloy 
thermodynamics is enforced through the Ginzburg-Landau functional 



(6) 



where the coefficients can be expressed as r' = ks(T c — T)N V and c = kBT c N v \^/2, with II the critical temperature 
of the binary alloy and Ao the interaction distance which is of the order of lattice sparingly In the absence of elastic 
effect, when T > T c the equilibrium state of bulk alloy is homogeneous with = 0, while for T < T c two phases 
4> = ±y/r' /u coexist in equilibrium. The gradient energy term in Eq. @ is important for stability analysis. It 
penalizes the-skarp compositional changes and without this term a nonphysical divergence for short wavelength mode 
will occurHtatL3 The coefficient c of this gradient energy term is related to the interfacial tension arising at the 
interface between coexisting phases. 

The elastic free energy functional J- c \ is important in any stress-driven system such as a growing film. Linear 
elasticity theory provides its expression as 
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For the 



where cr^ is the stress tensor and Syjw is the elastic compliance tensor with subscripts i, j, k, or I 
isotropic system assumed here, we have Sijki = Sik$jl0- + v)jE 
ratio v. Generally, the elastic constants of the film depend on the local composition. In the first order approximation 
they are 



SijSkiv 'I ' E, with Young's modulus E and Poisson 



H f {l + fi<j>) 



(8) 



and v* — E> lln* — 1 (with \i the shear modulus), and the composition-independent results can be obtained by 
setting E\ = fi* = 0. In this paper, superscript / or s refers to the film or substrate respectively, and in most of our 
calculations below (except Appendix we neglect the difference between the average elastic constants of the film 
{Eq, fiQ, and i>l = Eq/2(j,q — 1) and of the substrate {E s , fi s , and v s ). This corresponds to the material systems with 
substrate and film having similar elastic constants. 

Using the model developed here and the equations (Q)-(||), we can explore the stability of the growing symmetric 
alloy films in the parameter space of lattice mismatch e, solute expansion coefficient 77, deposition rate v, growth 
temperature T, and the composition dependent elastic moduli parameters E\ and fi\. 

III. SOLUTION OF MECHANICAL EQUILIBRIUM 

Thermodynamic equilibrium consists of mechanical, thermal and chemical equilibria. In the nonequilibrium system 
of growing alloy films that we are considering, the mechanical equilibrium occurs on a very fast time scale and we can 
assume that it occurs instantaneously. Thus by solving the equations of mechanical equilibrium, we can obtain the 
elastic displacement field u in terms of the other two fields <j> and h. This is the goal of this section. Knowing this 
solution, its substitution in the elastic free energy functional J- c \ renders it dependent only on <j) and h. The resulting 
effective elastic free energy functional is then used in the dynamic equations (J2J) and (^) for further analysis. 

Inside both the film and substrate, the mechanical equilibrium yields 



djdij = 0, 



(9) 



where the linear elastic stress tensor er^ obeys the Hooke's law, and for the isotropic system studied here it is given 
by 



4 = v 



(e + rjcj>)5i 



(10) 
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in the film, subjected to both the misfit strain e and compositional strain rj<p, as well as 
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v 



1 - 2v 



■ufAj +u 8 ij 



(11) 



in the substrate. Here Uy = (diUj + djUi)/2 is the linear elastic strain tensor, and the reference state we use is the 
unconstrained state of substrate bulk lattice. 

Considering the negligible external pressure on top surface of the film and the continuity of displacement and stress 
at the coherent film-substrate interface, we have the corresponding boundary conditions: 



aijiij=0 at z = h(x,y,t), (12) 



with rij the unit vector normal to the surface, and 

a ij n j = a ij n j an( l U i — U i at z = C{ x ill), (13) 

where rij is the unit normal to the film-substrate interface and oriented toward the film. Moreover, in the substrate 
the displacement and strain are expected to vanish far from the interface, i.e. 

ttf — > and ufj — > for z-> — oo. (14) 

To solve the above mechanical equilibrium equation, all the variables are expanded in Fourier series with a general 
form 

q 

where £ could represent different variables such as 0, h, Ui, and £, and £(q, z, t) denotes the small perturbation around 
the basic state solution £ (for height variable h or interface variable £, the perturbation should be expressed as h(q, t) 
or C(q), respectively). In the basic state, we have a uniform growing film with homogeneous composition = and 
a planar front moving at constant rate v, corresponding to film thickness h = vt. Due to the coherency of the film, 
the in-plane film lattice constant is equal to a s and then u{, = = 0. The Poisson relaxation in the z direction 

results in u{ — u[ z z, with u{ z — u = e(l + )/(l — ). The film is stressed in the lateral directions, leading to 
al x = af y = a = — and all the other strains and stresses are zero. In the substrate, we have uf — and 

U'j = erf 7 - = for i, j = x, y, z, corresponding to an unstrained basic state 



The expansion forms of both the mechanical equilibrium (^) and boundary conditions (12 )— (|T4|) can be obtained 
using Eq. (p^| ). Generally, due to composition dependence of film elastic constants, as depicted in Eq. (||), from Eqs. 
(^) and (|10[) one can obtain the extra terms proportional to coefficients E* and /i* in the mechanical equilibrium 
equation of the film, which is different from the composition independent case and more complicated (since the elastic 
constants also vary with position r according to Eq. (ph). Here we only consider the solution of mechanical equilibrium 
to zeroth order of elastic moduli and , which is applicable when E\ and [i\ are small enough, or and (i* are 
assumed to vary very slowly with space and the equation is considered to first order. Consequently, we obtain the 
mechanical equilibrium equation with a linear form as used beforeE3, but with different linearized boundary conditions 
due to the film-substrate interface roughness (. At the top surface of the film, i.e. z — h, from Eq. (|l^) we have 

Hz = -iqxvh, 



a 



iq y crh, 



H z = 0, (16) 
as before, while at interface z — ( = 0, we have different conditions due to Eq. (|l3|): 

-iq x a( + Hz = Kz, 

f f = a s 

yz "yzi 



-iq y a( + Hz = K 



Hz = Kz, (i7) 

and 

u{ = u\ for i = x,y, 

uC + u{ = u s z . (18) 
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Finally, for z — > — oo condition (n4) leads to 



and d z ul 



(19) 



The corresponding results of the solution for non-planar interface £ ^ and differing film-substrate elastic constants 
are detailed in Appendix |a|. In-the following studies for elastic free energy and dynamical equations, we assume a 
planar film-substrate interface,!^ i.e. £ = 0, as well as the equal average elastic constants between film and substrate, 
that is, Eq = E s — Eq, ^ = [i s = fio, an d Vq = v s = vq. Thus, substituting these conditions in the results of 
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Appendix IaL one can obtain the solution 
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iq y /q 
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1 - Vq 



V 



iq x W 
iq y W 
d z W 



for the film, where q 2 = q 2 + q 2 , and the coefficients cti and C are given by 

-qh 



iq x at x + iq y a y 



(2ii - /'„) - ,jln [ qnh 4- ( i±^j m (d z W /z i 



1 - Vq 



z—h 



(20) 



(21) 



qa z = e-i> 1 



(l - 2^ + f guft + (jz^J m (d»W / 
-(2(l-v )+ q h)(\±^)r,(d 2 W- 

\X — 1Sq J \ / z—h 



(22) 



1 - z/ 



Tj (q d z \V — ■ 



z—h 



(23) 



with q y a x — q x a y . Here W is introduced through the particular solution of the mechanical equilibrium equation, as 
shown in the last term of Eq. (|20|) , and is defined through 



(d 2 z - q 2 )W = <j>, 



(24) 



or the relation W 2 W = <j) in real space.E3 Consequently, the results for strain and stress tensors can be obtained 
through e.g., Eq. (|a|). 



IV. DYNAMIC EQUATIONS AND PERTURBATION GROWTH RATE 



To evaluate the dynamic equations (p|) and (^) and then determine the evolution behaviors of surface profile, we 
should first calculate the total free energy T . The first two contributions, surface and Ginzburg-Landau free energy, 
are given in Eqs. (||) and (||) explicitly, while the film elastic energy can be obtained from Eqs. (0) and jl0|): 



d 3 r 



1 . t f 2 ft- 



1 - 2vf 



° 2 / 

2 £ " £U " 



3e?70 + -rj 



V<t>u>u 



(25) 



with the Lame coefficient \f = 2fi^v^/(l — 2i>J). In Eq. (p5|), the first two terms represent the general form 
of the isotropic elastic energy density for single-component system, and by setting e = and renormalizing the 
Ginzburg-Landau free energy, which absorbs the </> 2 term, we can obtain the bulk elast ic energy functional with the 
effect of compositional strain (i.e. coupling term Tj^ul) that has been widely used.EJtJ Considering the composition 
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dependence of both Young's and shear moduli, we can write Eq. ( j25| ) to first order of E* and (i* as 

-h 



f 



1 



S 



l-2i/o 



3^* - 2(1 + j/qK 



2(1-2^) 



3 ^ 



1 



1 2 , 3£ x * - 2(1 + v )nl 



l-2i/ 



677 (/) - + 



3^* - 2(1 + z/q K 
l-2^o 



2(1 - 2u ) 



E{ - 



1 + 2u 2 



/2 ZE{ - 2(1 + , 



l-2u 



3 / 3£ 1 '-2(l + i / K 



(26) 



1-22/0 

with Ao = 2/Jo^o/(l — 2^o). The corresponding Fourier expansion to second order of perturbation is given in Appendix 

The strain tensor in elastic energy expression ( |26| ) can be determined from the solution of mechanical equilibrium 
shown in Eqs. (f20|)-(p3|) or Appendix |A|, and is in fact composed of two parts: 



(27) 



The first part is the homogeneous solution and can be calculated through the first two terms of the displacement 
solution (EM). The second part corresponds to the particular solution and is given by 



u par = i + i^^^.W. 



l-l/Q 



(28) 



Note that the homogeneous part uf° m is in fact not the function of system variable (j> or W, and the coefficients in 

this part, oti and C in Eqs. (f2lj)-(p3|), depend on (dzW);,^, (d^W)^^, and 4>\ z =h which are surface values (of fixed 
z for certain time t) evaluated at z = h due to the boundary condition at top surface. Therefore, when performing 
the functional differentiation of the fpee energy with respect to system variable 0, only the particular solution ( p8| ) of 
strain tensor needs to be considered. 

The dynamical equation for compositional perturbation at the surface can be determined by substituting energy 
formulae (H), (m and (E6n into Eq. |}|). When expanded to first order, it is 



84> s /dt : 



\-r' + cq 2 



- A^ s 



(29) 



with the elastic energy term (W 2 SJ 7 ^ l /S(p) q \ s obtained from the Fourier expansion using Eq. (^ 



V 2 ^ 



= -q 



E (3-vq) 2 £ (2(5 - 7v )Ej - (7 - 10u Q )(l + up)^) 
j-v + r, — 7~~vTFt — — ^ 



L (i-^o) 2 



(1 - i/ ) 2 (l - 2i/ ) 



Eq E (2Et-(l + 2v Q ) f i* 1 ) 

T) + r~rz ; £ 

l-v 



(1 - 2i^)(l - mb) 



fi «l 



E 



1 - i/ 



V 2 <^ s , which implies that [(&. 



(30) 

W 



where we have used the continuity condition of V 2 </> at surface: (V 2 0) 

since we ignore terms of order h 2 <ft s . 

For the evolution of surface height variable h, we can obtain the following linearized dynamical equation by com- 
bining Eq. (||) and the free energy functionals (||), (||) and (p(f): 



dh/dt = -T h q 2 ( iq 2 h + £ f 



(31) 



where £f\ s is the first order elastic free energy density evaluated at the surface, with the expression given in Eq. ( ]B3| ) 
of Appendix H. 
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Since what we are interested in is the stability property of this growing alloy system with respect to small pertur- 
bations, we focus on the film behavior of the early evolution time regime, where the morphological and compositional 
perturbations at the surface are assumed to evolve exponentially: 



n h t 



h = h {q)e^ 



and 



^ s = Mo)** 



(32) 



with f2/j and f2^ the growth rates of surface perturbations for h and 4>, respectively. Due to the coupling of dynamical 



equations (|2|) and (]3j) in general case, the joint stability or instability, corresponding to = ^<t> — is obtained. 
Considering the expression ([!]) for bulk composition field, we haveEl cj> = <fo exp(f2z/u), and then W = v 2 (j>/(fl 2 — q 2 v 2 ) 
from Eq. (E4). Consequently, by substituting this approximate form of W into the solution for film displacement 
vector w{ (|20)-(|23]), we can calculate formula (3^) as well as the linear elastic energy density £f\ s in term of two 
surface fields h and cf> s , as shown in Appendix p. 

For further analysis, it is convenient to introduce a characteristic length scale representing the typical size of 
interfaces between surface domains: 



as well as a diffusive time scale 



r 2 



-1/2 



(33) 



(34) 



where r = r> - 2E r) 2 /(l - v Q ) = fc s (T c ctt - T)N V is a renormalized coefficient with the effective critical temperature 
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(35) 



Note that T£ s used here is just the bulk coherent spinodal temperature derived by CahnQ for binary stressed alloys. 
Using the following transformations 
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(36) 



where we have set Th/T^ = 6 = la, as well as 

h* = h/l , 

4>* s = 4>B, 

we reduce the dynamical equations ( jBlj) and (|29| ) for surface perturbations to dimensionless forms: 

k 2 



dh*/dr = (e* 2 fc 3 - j*k 4 )h* - 



e Tj 



a — kVv \ 
a + kV ) 



2E{ ~ (1 + ^ 
2(1 - vq) 



(37) 



(38) 



and 



l/dr = 



1-Vq 



-k z \k 2 ±l 



(1 - 2v )e*rj* + (2E* - (1 + i/ )Mi> 
8^-5(1 + ^)^1 



2(1 -v 2 ) 



-e rj 



- V 



(1 - 2v Q )rf 2 + (2E{ - (1 + ^oK)eV 



k 3 V 
a + kV 



(39) 
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where for sign "±" , the top sign applies when the system is above the effective critical temperature T^ ff and the 
bottom sign is taken if T < T c eff . 

Due to the rescaling (|36|), the early time behavior of perturbations ( p2| ) is transformed to h* = Kq exp(<r^T) and 
cj>* = 4>q exp((T0r) with the nondimensional perturbation growth rates ah and Generally, the dynamical equations 
( |38| ) and ( |39"| ) for h* and cj>* couple with each other, resulting in a cubic characteristic equation for perturbation growth 
rate Oh~o^= a: 



(cr + 7*fc 4 -e* 2 k 3 ) 



2(1 - vl) 



— ((1 - 2, )^ 2 + {2E { - (1 + voMVrT) ^ + 



(1 - 2u )e*r)* + (2E* - (1 + vo)fi)e* 



V 



a-kVvo\ 2Et-{l + v )ii\ , 2 
e ?7 I r^r- I H — ; e 



(7 + fcV 



2(1 - i/q) 



0. 



(40) 



When the real part of a is larger than zero, the initial perturbations will grow with time, leading to the instability 
of the growing film; otherwise the surface perturbations will decay in time, making the system stable. Note that this 
characteristic equation (^p| ) is asymmetric with respect to the sign of rescaled misfit e*, leading to the asymmetry 
of stability properties between compressive._and tensile films as specified below, but it is symmetric for (e*,rf) — > 
(— e*,— ry*), similar to the previous results.BE3 Thus, in the following analysis we keep rj* > and consider different 
signs of e*. 

Setting El = /j,* — in above deductions we can obtain the corresponding results for system with composition 
independent elastic constants. The corresponding equation for h* is the same as that for h* in Ref. but the 
equations for </>* are quite different due to the fact that here we consider the total energy of the film and the 
surface-bulk coupling, while in Ref. [2(] only the surface energy is used (see also Appendix g). Consequently, the 
characteristic equations for a are also different, and our Eq. (^) redxups to a cubic, which_is_alHO different from 
that of the previous related work carried out by Guyer and VoorhccsBQL30 and Spencer et aZ.HoCJ due to different 
models and approximations. Compared with Eq. (f40|), this cubic equation shows the symmetry with respect to the 
sign of e* for systems with composition independent elastic moduli. 



V. STABILITY DIAGRAMS 



The typical dispersion relations of Re(<r), which is the real part of joint perturbation growth rate tr, as a function 
of rescaled wavenumber k are depicted in Fig. [l]for systems with composition dependent elastic moduli and T > T^ ff . 
The value of Re(cr) in Fig. [I] for each k is in fact defined as the largest real part of the three solutions of characteristic 
equation ( |40| ) , since that will be the dominant growth in case of instability. Note that for very large k values we always 
have Re(a) — > — oo, pv f' T1 when 77* is large, corresponding to short wavelength stabilization. This is different from 
the previous workcHlilo where nonphysical divergence for k — > 00 and large r\* occurs due to the absence of gradient 
term in free energy functional. If Re(cr) is smaller than zero for all k, e.g., the dashed curve in Fig. |l], the system is 
stable; otherwise, we can find a band of wavenumbers for which the surface profiles are unstable to fluctuations. In 
the following we investigate the stability properties and the effects of misfit strain, compositional stress, deposition 
rate, as well as the composition dependence of elastic constants more systematically, by determining the maximum of 
Re(cr) with respect to all disturbance modes k for different material parameters. 



A. Analytic results for special cases 



We first consider some simple cases where the misfit strain or solute stress is absent and then the analytic results 
can be available. The simplest case corresponds to e* — rf = 0, i.e. a lattice-matched film with constituents having 
the same size, for which dynamical equatkms (pq ) and ( |39| ) decouple, resulting in different perturbation growth rates 
for h* and </>*, as found in previous workEl 



07, = -7 k , 

(j (p = -k 4 Tk 2 - V. 



(41) 
(42) 
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Formula (^Tj) implies the stability of film morphology since ah is always negative. Eq. fl42] ) yields the stable composi- 
tional profile for T > T° s (top sign "— "), whereas for T < T^ s the competition between the deposition of homogeneous 
materials and surface decomposition leads to a critical deposition rate V c = 1/4. When V > V c the film is stable with 
respect to compositional fluctuations, otherwise phase separation occurs. 

When misfit strain exists (e* ^ 0) but solute stress is absent (rj* — 0), we have different results for systems with 
composition independent and composition dependent elastic moduli. When E* =p, fi* = 0, the .height variable and 
composition field still decouple, and we can recover the result of Asaro and Tillerfcl and Grinfelda for morphological 
instability, which is 



Oh 



*2, „3 



7 



(43) 



The stability of compositional perturbation is also governed by Eq. ( f42[ ) , in accord with the previous result.E2l However, 
when the clastic constants arc dependent on film composition, E* ^ and [i\ ^ make the dynamical equations for 
h* and 4>* coupled, and then the joint perturbation growth rate 07, = a^ = a obeys a quadratic equation 



a + aid + <2q = 0, 



(44) 



with coefficients 



01 = 7*fc 4 - e* fc J + k\k A ± 1) + V, 
a = {~f*k 4 -e* 2 k 3 )[k 2 {k 2 ±l) + V} 



[21% - (1 + ^oK] 2 *4, 5 
2(l + v )(l-v )* 6 



(45) 



Note that the solution form for quadratic (|44|) is a = I —at ± \/a\ — 4aoJ /2, and we have maximum Re(cr) > due 

to a\ > and a < obtained from fl45| ) when fc <C 1, for both T above and below T®*. Therefore, for the lattice 
mismatched film without solute stress, the morphological profile and compositional mode couple with each other if we 
consider the composition dependence of elastic moduli, and then the misfit strain can not only cause the morphological 
instability, but also the phase separation in spite of uniform deposition, which is different from the previous stability 
resultsM 

So far in this subsection, we have considered different cases all of which have rj* = 0. From Eq. (|35|), we note that 
for these cases T^ s is same as T c . 

The case for no misfit strain (e* = 0) but nonzero solute stress (77* 7^ 0) can determine the role of compositional 
strain on film stability (without the interplay with misfit). The same results of perturbation rate can be obtained 
from dynamical equations for composition dependent (Eqs. d38| ) and (|39|)) and composition independent elastic 
moduli. That is, the morphological and compositional modes decouple, with ah given by Eq. ([flf), resulting in stable 
morphology of this lattice-matched film. The property of compositional perturbation is more complicated, with a^ 
governed by a quadratic similar to (|i^): 



where the coefficients are 



a\ + aia-0 + a = 0, 



ai 



no 



k 2 (k 2 ±l) + (k + l)V, 



W 



'.1,1.2 l 1 , 1.2 ( 1 \ - 



fe 2 (r ± 1) - k A 



1 - V 



n 



v 



(46) 



The stability conditions for compositional profile can then be calculated analytically through the solutions of this 
quadratic equation for a^. For growth temperature T higher than bulk coherent spinodal T° , which corresponds to 
the case of most experiments, phase separation can occur only when 77* is larger than a critical value 77*, while for 
r/* < 77* the lattice-matched alloy film is stable against decomposition. Different condition applies to low temperature 
T < Tf, for which there is no critical value like rj* for the appearance of alloy decomposition, and larger deposition 
rate V can suppress the instability. In more detail, the stability condition for T > T° ff is given as 



1 



l-2i/ 



#2 , * 2 

v < Vc 
or 

V* 2 >vf and V>\ 



for any V 



1 - whtr 



(47) 
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For T < T^ s the compositional profile is stabilized with the condition of high deposition rate 

1 r „n2 



V > 



1 + {r)*/r)* c 



(48) 



which is fulfilled for any rf , with the value of 77* given in Eq. fl47|). The corresponding stability diagrams calculated 
from Eqs. @ and @ are plotted in Fig. |, with v = 1/4 for which ?7* 2 = 1.5. In both Fig. | (a) and Fig. | (b), 
the compositional profile is stable in the region of high deposition rate V > V c , with V c generally increasing with the 
value of 77* except for a special region of T > T° s . This special region corresponds to 77* < 77* as given in Eq. (|47j), 
and in this range of rf we have V c — 0, that is, the system is always stable. 

Therefore, our resutts-for this instability driven by compositional stress, in the absence of misfit, are different from 
those of earlier works,"!! 2 ! where the compositional instability is expected only below an effective temperature which 
increases with ?7 2 , but is higher than the coherent spinodal T° s used here and independent of deposition rate V. 
Our calculations exhibit more complicated behavior, that is, the solute stress enhances the compositional modulation 
only when exceeding a critical value (77*) for T > T^ ff , and this stress driven instability can be suppressed by large 
deposition rate, as illustrated in Fig. g as well as Eqs. ( [47| ) and (|48|). Note that both the increase form of V c and the 
value of 77* in Fig. |^ are different from that of Ref. due to different models and characteristic equations for a^. 



B. Stability diagrams in general case 



Generally, both misfit strain and solute stress are nonzero in strained alloy film growth, corresponding to e* 7^ 0, 
77* 7^ 0, and the coupling of dynamical equations for surface height and compositional variables. Thus, the growth of 
morphological disturbance will induce the decomposition and phase separation of alloy film, and vice versa, resulting 
in the simultaneous occurrence of stability or instability for morphological and compositional profiles. The joint 
perturbation growth rate is then determined by the cubic characteristic equation (f40|), which can be rewritten as 



a + a,2<T + a\o + ao = 0, 

with the three coefficients: 

a 2 = 7 *fc 4 - e* 2 k 3 + fc 2 (fc 2 ± 1 + /fee*??*) + (k + 1)V, 



(49) 
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Here the parameters «o and [3q ar c defined as 



(1 - 2^o)eV + «oe* 2 -2vq(1 - v Q )e*rf + a e* 2 



(50) 
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8^* - 5(1 + up)t4 
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(51) 



We can determine the stability properties with respect to different parameters such as e*, 77*, V, and T from the 
largest real-part root of above characteristic equations, which can be obtained from the dispersion relations similar 
to those shown in Fig. [j] In this and the following sections two sets of material parameters are used, corresponding 
to two typical types of stability properties. In the first set, we choose vq = 1/4, 7* = 5, 77* around 0.6, as well as 
Ei = —0.4 and /z* = —0.1 if considering the composition dependent elastic moduli, which qualitatively represents the 
alloy films of SiGe type. The second set corresponds to the systems of pseudo-binary III-V alloys such as InGaAs, with 
the chosen parameters v = 1/3, 7 * = 3.5, 77* around 1.1, and E{ = —0.25, = —0.5 if applied. These parameter 
values have been arrived at by using the rescaling formulae j36| ) as well as different types of material parameters for 
group IV and III-V components presented in Ref. || (e.g., the elastic moduli and the atomic size difference for alloy 
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components that determines 77). The parameters such as the elastic constants and lattice parameters are estimated 
from the average of the quantities for pure components (e.g., Si and Ge or InAs and GaAs), while E* and [i\ are 
adjustable but should be qualitatively consistent with the variation trend of E and v between pure alloy components. 

In Figs. g and g we plot the stability diagrams of dimensionless deposition rate V versus misfit strain e* for 
different values of 77* , including different cases of growth temperatures (above or below coherent spinodal T° ) and of 
composition dependence or independence of elastic moduli. When one compares Figs. g (a) and g (a), the stability 
diagrams are qualitatively similar. For E* = /j* = 0, the stability boundaries are symmetric with respect to sign of 
e* and qualitatively insensitive to changes in vq or 7*. The major qualitative difference between these two figures 
is seen for the case of composition dependent elastic constants, i.e. when E* ^ and /zf 7^ 0, as we compare Fig. 
I (b) with Fig. I (b) which are for T > T° s . These stability results are clearly asymmetric with respect to the 
sign of misfit e*. Similar asymmetry also occurs for T < T c cff . This asymmetry emerges from the characteristic 
equation (|4^), and is enhanced by larger solute stresses (i.e. larger values of 77*). However, the form of asymmetry is 
dependent on the choice of parameters. In Fig. g, the parameters of SiGe-type are used and the stability diagrams 
show that the instability region of a film under compressive strain (e* > 0) is larger than that under tensile strain 
(e* < 0), while in Fig. g with InGaAs-type parameters, we have opposite result that the compressive filna^havej 
larger stability region. This, asymmetry of stability has also been found in earlier works of Guyer and Voorheesaall3'til 
and Leonard and DesailiSEj (for the case of composition dependent Young modulus) when single surface mobility is 
considered. However, in these results only one kind of asymmetry is obtained, that is, tensile films are more stable 
than compressive films, while our calculations yields different kinds of asymmetry due to different parameters E\ and 
\i\ determined by composition dependent behavior of elastic moduli. Note that most experiments are carried out 
at temperature T > T° ff , with growth parameters corresponding to very small rescaled deposition rate V and for 
rescaled 77* (~ 0.6 for SiGe film and ~ 1.1 for InGaAs growth) that is not large. Thus, to compare with available 
experimental observations, we should study the bottom part and the lowest curve of diagrams in Figs. I (b) and g 
(b), which, however, have very weak asymmetry. In the next section, we will present the results for kinetic critical 
thickness on the scales corresponding to experimental parameters, and the asymmetry of (effective) stability will be 
shown more clearly, with more complex behavior for InGaAs-type films. 

A common feature of the diagrams in Figs. g and gis that the system can be completely stabilized by large enough 
deposition rate V, the value of which is determined by material parameters like e* and 77*. This stabilization effect at 
large deposition rate can be understood as follows. First, when the materials are deposited fast enough, the former 
surface layer is buried and then frozen (due to miniscule bulk mobility) before the surface morphological undulation 
and phase segregation have enough time to develop through surface diffusion. Second, as described in dynamical 
equation (|j|), the deposition flux is uniform and then have the tendency to drive the surface composition profile to 
the homogeneous phase of incident flux. Note that in Fig. g the scale of V in the stability diagrams ((c) and (d)) 
for T < T^ s is much larger than those ((a) and (b)) for T > T° ff , due to the fact that in usual coherent systems, 
spinodal decomposition is apt to occur below the temperature Tj? , and then faster depositiqapiajieeded to suppress 
it. Thisjsffect of stabilization at high deposition rates has ibeecL found by Guyer and Voorheesall3'til and Leonard and 
DesaitSEl, but is absent in recent work of Spencer et alx3i3 for the case of equal mobilities of alloy components. 
However, in the results of Guyer and Voorhees, the stress driven instabilities can be suppressed only for films under 
tensile strain (e* < 0), while for compressive layers this stabilization by high V is not possible and then systems 
are always unstable. This is different from our results where films of both e* > and e* < can be stabilized for 
symmetric or asymmetric cases, as shown in Figs. g and g. 

Compared to single-component film, one of the additional and crucial factors for alloy film stability is 77, reflect- 
ing the oca nposi tional stress caused by different atomic sizes of alloy components. From the model of Guyer and 
VoorheesaiJ'HEllitiSjfDund that the compositional stress can stabilize the film modulations. However, recent study 
of Spencer et aZ£3u3'Efl shows that when one considers single atomic mobility, the coupling of solute strain (77) and 
film-substrate misfit strain e always enhances the instability and makes the alloy film more unstable than the single- 
component film, and only for the case of unequal surface mobilities can the solute strain lead to a decrease of the film 
instability with a possibility of stabilizing the system through deposition. Here we have used one effective surface 
mobility to describe the dynamics of the alloy film, but have more complicated and different results for the role of 
solute stress (or 77) . To obtain the effect of 77 more clearly, without the interplay of composition dependence of elastic 
constants, we consider the case of E* = pb* = 0, corresponding to Figs. g (a), (c) and g (a). These three diagrams 
show that the joint morphological and compositional instabilities can be suppressed by large deposition rate, as dis- 
cussed above, resulting in less unstable film compared with single-component system where the Asaro-Tiller-Grinfeld 
morphological instability (as described in Eq. (|42|)) always occurs for any deposition rate. In this sense our result is in 
agreement with that of Guyer and Voorhees. However, in the presence of solute stress the diagrams in Figs. g and g 
exhibit the smaller stable regions for larger 77*, corresponding to the destabilizing effect of large compositional strain, 
which is opposite to the effect found by Guyer and Voorhees but similar to that of Spencer et at. A new feature of 
our results is in the stability diagrams for T > T£ s , i.e. Figs, g (a) and g (a), where the system is always stable for 
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small misfits |e*| < £q when r\* < 77*, with 77* defined in Eq. (E7[) and eg dependent on the value of 77*. Thus, only 
when the compositional stress exceeds a critical value, i.e. 77* > rj*, can the instability of the film be significantly 



enhanced. This is different from all the previous results, and can be explained by Eq. (47) in the limit of zero misfit, 
which shows analytically that a film is stable against any disturbance when 77* < 77* . This feature does not exist for 
T < T° s , which can also be expected from the stability condition ( pf8| ) for e* = 0. 

VI. KINETIC CRITICAL THICKNESS 

According to the stability calculations shown above, in general the morphological and compositional modes are 
coupled, that is, instability in one also implies instability in the other, and then the joint stability or instability is 
determined by diagrams illustrated in Figs. |^ and ||. However, in some experiments, although the chosen material 
parameters and experimental conditions correspond to the unstable region of our theoretical diagram, the film profiles 
are in fact obs er . ye d to be stable. This may be due to the kinetic stabilization effect of growing film, as proposed by 
Spencer et aLPir 6 ! When the real part of perturbation growth rate a is positive, the instability develops and grows 
with time. However, at early time of growth the surface perturbations have not enough time to develop and become 
observable compared to constant deposition of materials flux. Only after a sufficient time, that is, when the film 
is thick enough, can the perturbations develop substantially and be apparent relative to the planar basic state of 
the film. Then the perturbation growth rate (a) can exceed the relative growth rate of the film basic state (v/h). 
Thus, this competition between deposition and perturbation growth results in a kinetic critical thickness h c , with the 
nondimensional form given by 

K = V/a (52) 

in linear analysis Jlifcl 111 where a is usually approximated by the maximum perturbation growth rate cr max with 
respect to different growing mode k. Below this critical thickness /i c , the instability of the growing film is effectively 
suppressed by deposition flux and may not be seen experimentally, while above h c the surface undulations are seen. 
This phenomenon for the onsetpof stress driven, instab ility at certain film thickness has been observed by experiments, 
for both the systems of SiGeBEa and InGaAs.EE§EZl 

^From the mechanism of kinetic stabilization described above, one can expect that kinetic critical thickness h* 
will increase with deposition rate, which is verified by our calculations of h* versus V shown in Fig. for growth 
temperature T > T° ff . For both the SiGe-type and InGaAs-type parameters, the log-log plots in Figs. [5[ (a) and (b) 
show that when the rescaled deposition rate V is not large (corresponding to most experimental cases) , generally the 
relationship between h* and V is not linear, and the nature of the rise of h* with V for compressive and tensile misfits 
is different for the case of composition dependent elastic moduli. At sufficiently high deposition rate, the curves of 
h* for opposite signs of misfit e* and for different composition dependent forms of elastic moduli converge, and its 
variation with respect to deposition rate is linear, which can be seen on scale of V larger than that of Fig. ^|but is not 
shown here. These properties imply that the maximum of perturbation growth rate a depends on small deposition 
rate V as well as on the values of E*[ and but is independent of them when V is large enough. 

^From the plots in Fig. [| we can also obtain the results of stability asymmetry with respect to the sign of misfit, 
except for large enough deposition rate V which results in the presence of symmetry of h c in e. For SiGe-type system 
with composition dependent elastic moduli, Fig. |5| (a) shows that kinetic critical thickness for compressive films 
(e* > 0) is smaller than that of tensile films (e* < 0) when V is not large enough, corresponding to. the phenomenon 
that the growing films under tension are effectively more stable, as observed in SiGe experiments.^ The asymmetric 
property for InGaAs-type system is more complicated, that is, it depends on the deposition rate for T > T° ff , as 
shown in Fig. [5] (b). For intermediate values of V, the compressive films have larger h* compared to the films under 
tension, making compressive films more stable. For small V, h* of compressive films is smaller, leading to less effective 
stability which is similar to SiGe-type films. Note that the asymmetric property obtained here for InGaAs is different 
from the result of higher 77* values shown in the stability diagram of Fig. ^ (b), which in fact represents the asymmetry 
of absolute stability, but not that of effectively kinetic stabilization discussed here. 

This deposition dependent behavior of asymmetry for InGaAs films, can qualitatively explain the contradiction in 
recent experimental observations. In the experiment of Okada et al.f^i conjipressively strained layers are found to be 
more stable than tensile layers, whereas the observation of Guyer et al£3 yields opposite result. According to our 
calculations here, this inconsistency is attributed to different deposition rate used in these experiments: The growth 
conducted by Guyer et al. has much lower deposition rate and then corresponds to the asymmetry similar to SiGe 
films but opposite to that of Okata-ef al.. The direct measurement of critical thickness for the onset of instability has 
been carried out by Gendry et al&3 for different signs of misfit and different deposition rate. Since our calculations 



mainly focus on 50-50 mixture, different from the alloy concentrations used in Ref. 35, here we only compare the 



theoretical and experimental results qualitatively. As observed during MBE growth of Gendry et at, for the undoped 
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layers with As-stabilized surface showing strong (2 x 4) reconstruction, the critical thickness in compression is found 
to be less than that in tension when deposition rate is small, while for high deposition rate, larger h c in compression 
is obtained. This is consistent with our findings in Fig. || (b) for the case of composition dependent elastic moduli. 

The scaling behavior for critical thickness h c with respect to misfit e has also attracted attention in recent work. 
Since the system becomes more unstable with the increase of misfit strain, h c is expected to decrease with larger 
e. For single-component, lattice mismatched ifijw, where the effect of solute stress is absent, the kinetic critical 
thickness is found to follow the e~ 8 power lawJl§E3 Our calcuations also lead to e -8 power law as follows: From Eq. 
([52|), h* — V/a max ; from Eq. (|43|), one can show that the maximum growth rate occurs at fc max = 3e* 2 /Aj* which 
leads to h* — (2567* 3 V)/(27e* ). Th&-pbserved behavior for multi-component system like SiGe is different. The 
experiment carried out by Perovic et alEB finds that kinetic critical thickness h c for coherent SiGe/Si(100) varies with 
misfit as e -4 , while recent observation of Tromp et alr3 found a much slower dependence e^ 1 . Our results for the 
dependence of rescaled critical thickness h* on e* are shown in Fig. ^, and indicate that generally the scaling of h* 
versus e* is not a power law. For the SiGe- type system, as illustrated in Fig. || (a), the curve (dotted) for the case of 
composition independent elastic moduli {El = /i* = 0) presents a behavior close to 1/e* 8 for large misfit e*, but has 
slower decreasing form for small e* . When we consider the composition dependence of elastic constants {El ^ and 
//* 7^ 0), the decrease of h* with e* is slightly faster than e*~ 8 for the compressive layers (e* > 0), but much slower 
for the tensile films. In Fig. ^ (b), corresponding to the InGaAs-type case, h* scales similar to e*~ 8 only for small 
misfit e* , but decreases much slower at larger values of misfit, in particular for films grown under tension if El =/= 
and fx\ ^ 0. 

The asymmetry of effective stability {i.e. thickness h c ) with respect to the sign of misfit can also be seen in Fig. 
for the case of composition dependent elastic moduli. For the SiGe- type system (Fig. || (a)) h* for tensile films 
is always larger than that of compressive films, that is, tensile layers are effectively more stable, as found earlier. 
However, the InGaAs-type system yields more complicated scaling behavior. For large misfit e*, Fig. || (b) exhibits 
that films under tension have higher h*, and then can be observed as more stable, similar to the result of Fig. |6| 
(a). This is in accord with the result of h* versus V shown in Fig. || (b), since here the rescaled deposition rate V 
is chosen as 10 -3 , corresponding to the small V region of Fig. |] (b) where one has large |e*| = 0.8 and SiGe-like 
behavior of asymmetry. For intermediate values of |e*|, asymmetry opposite to SiGe-type system is obtained, that is, 
the h* curve for compressive layers lies above that for tensile layers, resulting in more stable compressive films. When 
|e*| gets smaller (not shown here), the difference of h* values between e* > and e* < becomes negligible, and 
then we have the symmetry of stability for films under compression and under tension even when elastic constants 
depend on composition field. Therefore, combining Figs. || (b) and|| (b), we conclude that the asymmetry of effective 
stabilization for InGaAs-type films is dependent on both the deposition rate and the magnitude of misfit strain. 

VII. CONCLUSIONS 

In this paper, we have determined the stability properties for growth of dislocation-free and coherent, strained alloy 
films by developing a continuum dynamical model and performing a linear stability analysis. In our model there 
are two important and new features, the coupling between top surface of the film and the bulk underneath, and the 
composition dependence of elastic moduli {E and fx), which result in new stability results. Both the thermodynamics 
and elastic effects are considered. Besides the phase separation mechanism at low temperature and the role of surface 
energy that favours the planar surface, we also include the effects of elasticity which can arise due to (i) the mismatch 
of lattice constant between a growing film and the substrate, resulting in nonzero misfit strain e; (ii) the dependence of 
lattice constant on composition, leading to nonzero solute strain represented by a coefficient 77; and (iii) the dependence 
of elastic moduli on composition, as given in Eq. (Q). Moreover, the deposition rate can play an important role for the 
growing films. The main results, due to the interplay of these factors, that we obtain can be summarized as follows. 

For single-component epitaxial growth, when one studies the stability of the free surface of a growing film, the 
morphological corrugation may lead to the formation of coherent islands. This occurs on account of the misfit strain 
caused by lattice mismatch. In contrast, for multi-component growth like MBE, the time evolution of the morpho- 
logical profile is coupled to the dynamics of composition. Physically the coupling occurs since smaller (larger) atoms 
of the film prefer to incorporate in regions of compressive (tensile) stress. This leads smaller atoms to preferentially 
diffuse to troughs (peaks) and larger atoms to segregate to peaks (troughs), for a film under compression (tension). As 
a result of the coupling, there is a joint modulation of the surface morphological profile and the alloy composition, as 
well as a common perturbation growth rate a (see Eqs. (|38|), (|39|), and (ffOj)). Thua.a planar surface /ran be stabilized 
only if alloy decomposition is suppressed and vice versa. Some experiments (SiGeo and InGaAsEiJHo) have verified 
such a simultaneous development of morphological and compositional modulations. 

For bulk strained alloys, the compositional profile is homogeneous above the bulk coherent spinodal temperature Tj? s , 
and the decomposition can occur only below T° . However, for multi-component epitaxial films, alloy decomposition 
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occurs both theoretically and experimentally even for temperature T higher than T° s . This is due to two major 
effects. First, due to the coupling between morphological and compositional profiles, the phase separation can be 
driven by the difference in strain energy densities along the surface caused by the morphological undulation. Even in 
the absence of 77, that is, when the atomic size difference is negligible, the composition dependence of the film clastic 
constants can be the cause of coupling between morphology h and composition 0, as shown in Eqs. ( fli] ) and p3) for 
special case of e* 7^ but 77* = 0. The second effect is related to solute stress (or 77) itself. As depicted in Eqs. (|46| ) 
and ( |T"i| ) for e* — and 77* 7^ 0, when the morphological and compositional degrees of freedom are decoupled and we 
have stable surface morphology, alloy decomposition can also occur for T > Tj? if 77* is large enough. This means that 
when the solute stress is too strong, corresponding to very large difference in the atomic sizes of alloy components, 
the small and large species prefer to segregate with respect to very small background disturbances of surface profile. 

This compositional instability, as well as the morphological modulation, can be completely suppressed by high 
deposition rate, both for T above and below coherent spinodal T£ s , as seen in the stability diagrams of Figs, and 



|4| as well as from the analytic results derived for special cases in Sec. VA. Generally, increasing the magnitude 
of misfit e increases the instability, and larger solute strain (77*) also enhances the instability for most of the cases. 
Furthermore, the nonphysical short wavelength instability is absent in our work, due to the consideration of gradient 
free energy. 

Including the coupling between film surface and underlying bulk has resulted in some interesting results. One of 
the new features is on the role of 77. Our analysis naturally brings a critical value 77* for the stability at temperature 
T > Tf; . For 1 77* I below 77* , the system is always stable within a range of misfit e*, regardless of any deposition rate 
V. Only for \rj*\ > 77* can the increase of 77* significantly enhance the instability; this is shown in Eq. (4?]) and Fig. 
H for e* = 0, as well as in Figs. I and I for nonzero e*. When the surface-bulk coupling is neglected, the stability 
results do not have any critical 77*. Our results for the scaling behavior of kinetic critical thickness h c (the growing 
film is effectively stabilized by deposition up to thickness h c ) are also different from earlier work. In general, the 
dependence of nondimensional critical thickness h* on rescaled deposition rate V is not linear for intermediate and 
small V (corresponding to most experimental cases), and its dependence on rescaled misfit e* is not necessarily a 
power law, as shown in Figs. ^ and ^ respectively, even for = [i\ = 0. 

The asymmetry of stability with respect to the sign of misfit e* is found in the systems with composition dependent 
film elastic moduli, while for the composition independent case (E* = fi* = 0), it is absent. Figures || and || show 
the diagrams for absolute stability. For SiGe-type system, the films grown under tension (e* < 0) are expected to 
be more stable than those under compression (e* > 0), while the form of asymmetry is opposite for InGaAs-type 
films. For the effective stabilization, as measured by critical thickness h c , the asymmetry can also be found unless the 
deposition rate is too large. According to our calculations, a tensile film, with SiGe-type parameters, will be observed 
as stable up to a thickness larger than that for a compressive film, as shown in Figs. ^| (a) and || (a). However, for 
InGaAs-type system the asymmetry depends on both the deposition rate and the magnitude of misfit. As found in 
Fig. H (b), the system will have asymmetry similar to SiGe-type at small V, but for intermediate values of V reverse 
asymmetry is obtained, and then the compressivje-films. are effectively more stable. This can explain the observations 
of recent experiments on InGaAs epitaxial films P 1 tP 4 H 35 l When the deposition rate is fixed, the forms of asymmetry is 
opposite for large and intermediate magnitude of misfit e*, but at small misfit the symmetry with respect to the sign 
of e* is obtained (see Fig. || (b)). 

As discussed above and shown in Figs. 0-0, the stability properties and asymmetry depend on the chosen material 
system and the experimental conditions (such as deposition rate v and growth temperature T). E.g., if we desire a 
system with very strong modulations, stability diagrams (Figs. || and ^) and the variation of kinetic critical thickness 
(Figs. and |6|) show that the instability requires large misfit e, large solute coefficient 77, small deposition rate v, 
and large enough film thickness. Also note that the stability results presented here are with respect to dimensionless 
parameters V, e*, 77*, and 7*, which are all temperature dependent according to Eqs. (|36|), (^3|), and (|34j), via the 
coefficients r and T^,. The explicit temperature dependence of = T^/S is given after Eq. (0). Thus, even when 
we already choose the growing system with fixed material parameters, the stability behavior will vary with growth 
temperature, through the variation of the dimensionless parameters. 

All of the above results show that the problem of multi-component epitaxial growth, which is essentially a problem 
of nonequilibrium, is rich in interesting physics. It involves competition between thermodynamics and deposition, and 
coupling between morphology and composition. More interesting and richer results can be obtained if one further 
studies the detailed patterns and structures of the growing film, beyond the early evolution regime considered in this 
paper. 
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APPENDIX A: GENERAL RESULTS FOR MECHANICAL EQUILIBRIUM EQUATION 



constants, and then it is linear with the expansion form: 



As discussed in Sec. II], the mechanical equilibrium equation for the film is considered to zeroth order of elastic 
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with W defined by Eq. (|2J) of Sec. [11] and 
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The equation for the substrate displacement is similar to Eq. ( Al) but with <f> equal to zero. Considering the boundary 

(A4) 



condition ([L9|) for the decay of strains far from the interface, we have the solution form 
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which are used in the elastic free energy calculation of Sec. IV 
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for the substrate solution. In equations ( A7)-( A13 ), 
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To check whether these results are correct, we compare them with the known solutions by taking appropriate limits. 
For the single-componenLsystem with planar film-substrate interface and different elastic constants between film and 
substrate, Spencer et al.Ej have given the corresponding solution, which can be recovered by setting £ — as well as 



all the terms related to 4> and W to be zero in the above formulae (A7)-(A16). While for binary alloy system {(j) ^ 0) 
but with identical elastic constants for film and substrate as well as flat interface, the solution can be obtained by 
setting p — 1, 4 — vS — v 0i e ^ — e l; ^/ = b s , cf — c s , pf = p s , and C = in the above equations. The corresponding 
results are presented in Eqs. (pfj|)-(23) of Sec. Ill, and consistent with that of previous work.Ej 



APPENDIX B: ELASTIC FREE ENERGY IN FOURIER EXPANSION 

Using the Fourier expansion form (|l^), the elastic free energy functional ( p6[ ) can be expanded as 
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which has no contribution to dynamic equations (||) and (0), and the first order term 
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which is used in evaluating the linear form of Eq. (j^) for the surface height evolution. For the purpose of performing 
the linear analysis on Eq. (||) for the evolution of tj>, we should use the second order elastic free energy T^: 
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which is different from the result in Ref. ^ where surface-bulk coupling of the film is neglected and the second order 
clastic energy is evaluated at top surface z = h. If we calculate the functional differentiation ST^/Stp in Fourier space 
(that is, ST^I '<50(— q, z, t)) using Eq. (B4), we can also obtain the formula (^0|) of Sec. IV derived through real space 
procedure. 



The expressions for first order elastic energy density ( |B3| ) evaluated at the surface and the functional differentiation 
Eq. ( |30| ) in terms of h and <f> s can be obtained by using the strain tensor formula ( A6 ) of Appendix |A| and substituting 
the form of W derived in Sec. IV, i.e. W = v 2 (f)/(£l 2 — q 2 v 2 ), into the solution of mechanical equilibrium equation 

j L = 4>s to first order). They are given by 



and 



3|) (with the approximation . 

E Q 



I-Vq 



-2(1 + u )eqh + [2r) 



is qv 



2E* 



(1 + vq)h\ 



qv 



1-Uq 



(B5) 



,5T, 



cl 



= q 



2E r n + u 

l-vo \ \ 1 - v 



q [(1 - 2v )er) + [2E{ - (1 + ^oK)e 2 ] h 



a , 8Ef - 5(1 + i/ ) M i 
r? H — r- r e?7 



1 - VQ 



2(1 - v Q ) 
[(1 - 2*, )r, 2 



(2^-(l + ^K)e»7] 



qv 



9. 



qv 



(B6) 



respectively. Note that for the case of composition independent elastic moduli (E* = fl* = 0), Eq. ( JB5| ) yields a 
result for first order elastic energy density which is the same as that in Ref. The differences occur in Eq. (B6) 
due to differences in as well as different model and derivation procedure: In the previous work of Ref. 20, the 
second order elastic energy (B4) is evaluated at surface z = h first and expressed in terms of h and ^>| z= ft, and then 
the functional differentiation is carried out directly on surface <p with an approximation that 5Tl x /5cf> s ~ SF^/Scj)]^^. 
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FIG. 1: Dispersion relation for real part of perturbation growth rate a versus wavenumber k. 
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FIG. 2: Stability diagrams of rescaled deposition rate V versus 77*, for special case of no misfit e* = 0, vq — 1/4, as well as 
growth temperature (a) T > T^ s and (b) T < T° ff . These diagrams are calculated from the analytic results given in Eqs. (Ejj) 
and (El. 
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FIG. 3: Stability diagrams of deposition rate V versus misfit e* , with SiGe-like parameters v$ = 1/4, 7* = 5, and different 
cases of composition independent ((a) and (c): E* = n\ — 0) or composition dependent ((b) and (d): E{ = —0.4, nl — —0.1) 
elastic moduli, (a) and (b) correspond to T > T° ff systems, while (c) and (d) are for T < Tf. Stable and unstable regions 
are marked as "S" and "U", respectively. 77* values increase from bottom to top stability boundary lines, as indicated in the 
figures, and in figure (a) all curves with 77* < rj* = 1.5 1 / 2 cannot be distinguished from each other for the scale plotted here. 
Note that the vertical scales of V in (c) and (d) for T < Tf 1 are much larger than those of (a) and (b) for T > T^ a . 
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FIG. 4: Stability diagrams of deposition rate V versus misfit e* similar to Fig. except for T > and for InGaAs-like 
parameters vq = 1/3, 7* = 3.5, 77* = 2 1,/2 , as well as E* = —0.25 and /ij = —0.5. Compared with Fig. H, different asymmetry 
with respect to sign of misfit e* is shown for (b), due to different material parameters. The diagrams for T < T^ s (not shown 
here) are similar to those of Fig. til except the opposite asymmetry for nonzero E* and fi* . 




FIG. 5: Kinetic critical thickness (nondimensional) h* as a function of rescaled deposition rate V for growth temperature above 
Tf. The solid (e* > 0) and dashed (e* < 0) curves are for the case of composition dependent elastic moduli, while the dotted 
curve represents the result for composition independent moduli. The dot-dashed line of linear relationship h* oc V is also 
drawn for comparison, (a) The material parameters are chosen to qualitatively represent the SiGe-like film: uq = 1/4, 77* = 0.6, 
7* = 5, E* = —0.4, fi* — —0.1, and |e*| = 0.5; (b) The parameters are expected to qualitatively represent the InGaAs-like 
alloy: v = 1/3, rf = 1.1, 7* = 3.5, E{ = -0.25, fi\ = -0.5, and |e*| = 0.8. 




FIG. 6: Rescaled kinetic critical thickness h* as a function of absolute value of rescaled misfit strain |e*| for temperature 
T > . Similar to Fig. the solid and dashed curves are for positive and negative misfit e* respectively, with composition 
dependent elastic moduli, while the dotted curve corresponds to the case of composition independent moduli. The dot-dashed 
line represents the power law scaling behavior with exponent —8, which is used for comparison. In (a), the rescaled deposition 
rate is chosen as V — 10 -6 , and in (b) it is V = 10 -3 . All the other paramters are the same as those in Fig. 0. 



